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SUMMARY 


A new algorithm, which will transfer an aircraft from a given initial 
position and heading to some final position and heading, has been developed 
for the on-board synthesis of horizontal flightpaths. These so-called capture 
trajectories have the basic form of a turn along a circular arc followed by a 
segment of straight flight and a subsequent circular turn. The two circular 
arcs may have different radii. The algorithm finds all such solutions possi- 
ble, of which at least two and at most four exist, and selects the one with 
the minimum path length . Degenerate conditions in which one or more of the 
basic segments is missing are handled without difficulty. The solution to 
this problem is derived, and a Fortran listing of the algorithm is provided. 


INTRODUCTION 


Ames Research Center has been conducting studies involving use of on- 
board synthesis of flightpaths for area-navigation-equipped aircraft operating 
in the terminal area. The first step in the synthesis of a complete trajectory 
is the generation of the ground track. A useful criterion for generating the 
paths is to minimize the ground-track distance. The minimum-distance ground 
track, if flown at constant ground speed, has been shown by Erzberger and Lee 
(ref. 1) to consist of a series of straight lines and circular arcs. 

Pecsvaradi (ref. 2) developed algorithms for the on-board synthesis of hori- 
zontal flightpaths in that form. The flightpaths (e.g., as shown in fig. 1) 
consist of a fixed portion estaDllshed by predetermined waypoints and a 
capture trajectory that will transfer the aircraft from some arbitrary initial 
position and heading to the desired final position and heading defined by the 
capture waypoint. The basic form of the capture path is a turn followed by a 
segment of straight flight and subsequently a final turn. 

Flight- test experience has shown that the horizontal capture paths 
derived from the minimum distance criterion provide efficient and operation- 
ally acceptable terminal area trajectories for a variety of applications 
(refs. 3 and 4). The capture feature has been the subject of many favorable 
comments by pilots; however, it has been found that when the aircraft and 
capture waypoints are separated by less than four turning radii, the algorithm 
of reference 2, used in earlier experiments, often failed to find a capture 
path even though it is known (ref. 1) that at least two such paths always 
exist. 

The new capture algorithm developed here always yields the minimum- 
distance horizontal path of the basic turn-straight-turn form, including 


degenerate cases in which one or more segments of the basic capture path are 
missing. The set of possible solutions for the minimum-distance capture path 
(ref. 2) is not restricted to the basic turn-straight-turn pattern but may 
also include a sequence of three circular arcs. The latter option applies to 
a fairly restricted region in which the initial and final points are close to 
each other; it is eliminated from the algorithm in the interest of simplicity 
(as was done in ref. 2). The following section describes the algorithm and 
presents the derivations of the equations used in it. The additional equa- 
tions needed for implementation of the three-turn case are derived in appendix 
A, and a Fortran listing of the algorithm is given in appendix B. 


ANALYSIS 


Figure 2 is used to explain the problem and to define the variables. The 
turns are arcs of the circles shown in the figure and the straight portion of 
the trajectory must be a line tangent to both circles. Since the Initial and 
final turns may be either clockwise or counterclockwise there are four possible 
combinations of turning directions, two with the initial and final turns in 
the same direction and two with the turns in opposite directions. Figure 3 
illustrates one solution of each type. If a given pair of circles is entirely 
separate — that is, if no part of one circle lies within the other — it is 
possible to draw four tangent lines between the pair, but for only one of the 
four will the vector T) along the tangent line from the initial to the final 
circle coincide with the direction of rotation at both tangent points, as 
shown in the figure. Thus, there are at most four real solutions of the turn- 
straight-turn form. The conditions for a real solution are derived next with 
the aid of figure 2. 

For convenience, the vector TJ from the center of the initial turn to the 
center of the final turn is horizontal in the figure; however, the coordinate 
system is arbitrary. Furthermore, the location of the centers of the circles 
relative to an earth-fixed coordinate system— and hence the direction and 
magnitude of —depends on the direction of the turns. The constants and 
variables have been chosen so that the following discourse applies to all 
possible combinations of turning direction. 

Figure 2(a) is for the case where both turns are in the same direction 
and the tangent vector U does not cross in figure 2(b) the turns are in 
opposite directions and U crosses TJ. Initially the aircraft is at (Xi,Yi) 
in some inertial Cartesian coordinate system with heading Hj, defined as 
positive clockwise from the X-axis, and Vj is a unit vector in the direction 
of the velocity. The vector distance from (Xi,Yj) to the center of the turn 
is given by XT; Rj , where R; is the radius of turn and uj is a unit vector 
normal to V* and positive to the right of vi . Therefore, the vector from 
(Xj,Yi) to the center (XC lf YCi) is R;ui for a right turn and -Riuj for a 
left turn. The directions of turn are accounted for by writing the radius 
vector as RiSjui where S; ■ + 1.0 for right turns and Si « - 1.0 for 
left turns. Similarly the direction of the final turn is denoted by S 2 . 
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The aircraft moves along the circle from (Xi,Yi) to the tangent point 
(*2*^2) » which has a radius vector R^Si^. Vector If is the tangent vector 
from (X2.V2) at the end of the initial turn to (X^.Ya) at the beginning of the 
final turn. The radius vector at (X3,Y3) is R2S2U;* but since U2 and U3 
must be normal to D, U2 ■ U3. Likewise, the headings, H2 and H3, at the two 
tangent points are equal. The final turn ends at (X^.Yt,) with heading H4 
and radius vector R2S2U4. 

Using this notation we can write 

B + R2U2S2 “ R1U2S1 + TJ 


or 


^ ■ 5 + U2 (R2S2 " RiSj) ( 1 ) 

and therefore, since 1) and U2 are perpendicular, 

D - [Q 2 - (R*S 2 - RiSj ) 2 ] 172 ( 2 ) 

where by definition 

Q - [(XC 2 - XCj) 2 + (YC 2 - YC 1 ) 2 ] 1/2 (3) 

It can be seen from equation (2) that no real solution exists if 
Q < IR2S2 - RiSi|. When the turns are in opposite directions Si ■ -S2, and 
there is no real solution for Q < (Ri + R2)* that is, if the circles inter- 
sect. On the other hand, for rotations in the same direction Si • S2, and a 
real solution exists unless Q < | R2 ~ Ri I » that is, unless one circle lies 
entirely within the other. 


If a real solution exists, then from the definition of the radius vectors 

( — Rl S 1 sin HA 

I (4 > 

RjSj cos H\/ 


and 



/X Ci 

-Xi' 

\ 

RlUjSi • 





\YCi 


1 

Equating equations (4) and (5) gives 




XCi - Xi - 

RlSi 

sin 

Hi 

YCi - Yj + 

RlSi 

cos 

Hi 


(5) 


( 6 ) 
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Similarly 


XC 2 ■ Xif * R 2 S 2 sin H4 J 

(7) 

YC 2 ■ Y4 + R 2 S 2 cos H4 1 


The radius vectors at the tangent points can be used 
compute the components of Y 2 and 3f 3 

in the same manner 

X 2 - XCi + R1S1 sin H 2 

(8a) 

Y 2 - YCi - R1S1 cos H 2 

(8b) 

X 3 - XC 2 + R 2 S 2 sin H 2 

(9a) 

Y 3 - YC 2 - R 2 S 2 cos H 2 

(9b) 


Subtracting (8a) from (9a) and (8b) from (9b) gives the components of B: 
X 3 - X2 - XC 2 * XCi + (R 2 S 2 - RiSj) sin H 2 | 

Y 3 - Y 2 - YC2 - YCl - (R2S2 - R1S1) cos H 2 ( 

Another expression for the components of U is: 

X 3 - X 2 - D cos H 2 I 
Y 3 - Y 2 - D sin H 2 I 

Equating (10) and (11) gives 

D cos H 2 - (XC 2 - XCi) + (R 2 S 2 - RiS^) sin H 2 | 

D sin H 2 ■ (YC 2 - YCi) - (R 2 S 2 - R1S1) cos H 2 ) 

Equations (12) can be solved for 

(YC 2 - YCi)D - (R 2 S 2 - R1S1) (XC 2 - XCi) 
sin H 2 73 


( 10 ) 


(ID 


( 12 ) 


(13 


(XC 2 - XCi)D + (R 2 S 2 - R1S1) (YC 2 - YCi) 
cos H 2 - q5 


tan H 2 » 


(YC 2 - YCi)D - (R 2 S 2 - R1S1) (XC 2 - XCi) 
(XC 2 - XCi)D + (R 2 S 2 - R1S1) (YC 2 - YCi) 


(14) 

(15) 


Note that the quantities (XC 2 - YCi)/Q and (YC 2 - YCj)/Q are the sine and 
cosine, respectively, of the heading angle of 
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Equations ( 6 )- (9) and (13)- (15) completely specify a capture trajectory 
for one combination of Si and S 2 ; however, the length of the trajectory is 
needed in order to determine which of the feasible trajectories gives the 
minimum distance. The first turn angle is 


where 


The second turn is 


where 


TRi - (H 2 - Hi) + 2 ttCiSi 

1 0 if Si (H 2 - Hi) > 0 

1 if S 2 (H 2 - Hi) < 0 

TR 2 m (H 4 — H 2 ) ^ 2xC 2 S 2 

i o if s 2 (h 4 - h 2 ) i 0 

1 if S 2 (H 4 - H 2 ) < 0 


The lengths TD^ of the two arcs are 


TDi - Ri|TRi| 

td 2 - r 2 |tr 2 | 


(16) 


(17) 


(18) 


and the length D of the straight segment is given by equation (2), and the 
total length of the capture path is 

DT - D + TDi + TD 2 (19) 

It is possible to construct a "switching diagram" from which the flight- 
path yielding the minimum distance solution for any set of initial and final 
conditions is determined without finding all possible solutions. However, 
the switching diagrams are very complex, and from the computation standpoint 
it is simpler to solve for the length of all possible trajectories, as is 
done in the Fortran subroutine (appendix B) . 


APPLICATIONS AND EXAMPLE TRAJECTORIES 


The algorithm, which was Incorporated into an on-board guidance system 
for the Augmentor Wing Jet STOL Research Aircraft (AWJSRA), has been tested 
in simulations and during flight tests (ref. 4). In this application the 
fixed portion of the flightpath referred to in the introduction may Include a 
number of predetermined waypoints (e.g., as shown in fig. 3). The subroutine 
is called sequentially to synthesize the path between successive pairs of way- 
points working backward from the final one. The path between waypoints 3 and 
4 is synthesized by entering the following input variables 
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X! - X(VP 3 ) 

X* - X(WP 4 ) 

Yi - Y(WP 3 ) 

Y 4 - Y(WP 4 ) 

Rl - 0 

R 2 - 0 

Hi ■ 0 

H 4 - 0 

Note that if Ri » 0 then points (Xi,Yi), (XCi.YCi), and (X2,Y2> are 
Identical and the value of Hi is arbitrary; the same is true of (X 4 ,Y 4 ), 
(XC2,YC2>. and (X 3 ,Y 3 ) and the heading Hi*. The resulting trajectory moves 
from (X1.Y1) to (Xt,,Y 4 ) with heading 

1 Y4 

H2 * tan -1 rr- 
X4 

- Yi 

- Xi 

which (if waypoints 3 and 4 are specified properly) is the runway heading. 
The path from waypoint 2 to waypoint 3 is synthesized by calling the sub- 
routine with the input variables 

Xi - X(WP 2 ) 

Xi, - X(WP 3 ) 

Yi - Y(WP 2 ) 

y 4 - Y(WP 3 ) 

Hi - 0 

h 4 - H(WP 3 ) 

Rl - 0 

R2 • R2 (an input constant) 

As in the previous case, (Xi.Yi), (XCi,YCi), and (X2,Y2> are Identical and 
H4 is the heading at waypoint 3 computed in the previous step. The heading 
of the tangent path and the coordinates of the beginning of the turn are 


returned in H2 and (X3,Y3>, respectively. The synthesis of the fixed path 
is carried out in this fashion for successive pairs of waypoints until the 
initial waypoint is reached. Then the algorithm is used to synthesize the 
capture trajectory to one of the fixed waypoints. 

The photographs of the simulator cockpit display in figure 4 illustrate 
some close-in capture trajectories which would not have been found with the 
previous algorithm. If necessary for compatibility with ATC procedures, the 
algorithm could easily be modified to allow the pilot to select any of the 
possible capture trajectory patterns instead of the one giving the minimum 
path length. 


CONCLUDING REMARKS 


The solution of the tum-straight-turn capture flightpath problem 
developed here Involves generating all possible solutions, of which there are 
at least two and at most four, and then selecting the one with the shortest 
path length. Because all possible solutions are obtained, criteria other than 


m *■ 




minimum path length could ba uaad. For example, a nonminimum path length 
solution may ba chosen in ordar to avoid a raatrictad airspace region. The 
closed- fora aquations determining the solutions are easily solved on currently 
available area-navigation computers or on a microprocessor. Other Important 
flightpath problems in the terminal area can ba approached via the capture 
solution. These Include path stretching* holding patterns, and flight through 
a sequence of waypoints* all of which can be formulated as a sequence of turn- 
straight-turn capture problems. 
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APPENDIX A 


DERIVATION OF THE THREE CIRCULAR ARC PATTERNS 


Ic can be shown by geometrical construction that there are at most four 
turn- turn- turn type trajectories for the horizontal capture problem, namely, 
two each of the right-left-right and left-right-left patterns. However, it 
was shown in reference 1 that the middle turn must exceed ir radians for a 
minimum-path- length trajectory. This requirement eliminates one trajectory of 
each pattern leaving at most two eligible three-turn trajectories. 


The problem is illustrated in figure 5 for the right-left-right pattern. 
The vectors in the figure are defined aa before, recognizing that in this case 
S2 ■ Si. In order to satisfy the requirement that the middle turn exceed ir 
radians its center must be on the opposite side of 5 from the straight seg- 
ment of the rlght-stralght-rlght solution shown for comparison. Furthermore, 
no three-arc solution exists for Q > Ri + R2 + 2R3. 

Let Hq be the heading of 1$, defined in the text and define a unit 
vector vq with heading angle, Hq as follows 

H 0 - H q + Si | 


Then vq Is perpendicular to (J and points in the direction of flight where 
7 } intersects the circle of the initial turn. From the law of cosines 


cos A’ 


Q 2 + (Ri + R3) 2 - (R 2 + R 3 ) 2 
2Q (Ri + R 3 ) 


COS B' 


Q 2 + (R 2 + R3) 2 - (Rl + Rs) 2 
2Q (Ri + R 3 ) 


(Al) 

(A2) 


The direction of turn is accounted for by defining 

A - S^’ (A3) 


and 


3 - SiB* 

Using these definitions it can be seen from the figure that 

H 5 - H 0 + A 1 
He - H 0 - B + it I 


(A4) 


(A5) 


From the definition of Hq and the equations derived in the text 
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I r 


Equat: 

exprei 

Chang: 


■in Ho 


S!(YC 2 - YCi) 


coa Ho ■ 


-Si(XC 2 - XCi) 


(A6) 


(A7) 


(A 5 ) can be used with double-angle trigonometric identities to obtain 
tons for sin H5, cos H$, sin Hg and cos Hg, in terms of ll 0 . A, and B. 


(A8) 


appropriate subscripts 

in equations 

1 (8) and (9) in the text gives 

Xs 

- XCi 

+ R1S1 

sin H5 | 

Ys 

- YCi 

- R1S1 

cos H5 ) 

x 6 

■ XC 2 

+ R 2 Si 

sin % 1 


- yc 2 

- R2S1 

cos He ) 


The turn angles are calculated as follows 

TRi - (H S - Hi) + 2 nC 5 Si 


where 


C 5 


jO if (H 5 - H j ) S 1 ^ 0 
'l if (H S - H 1 ) S 1 < 0 


(A9) 


(A10) 


where 


TR 2 - (H4 - H 6 ) + 2 nC 6 Si 


C 6 " 


to if Oh* - h 6 )s > 0 
(1 if (Hi* - H 6 )S < 0 


(All) 


where 


TR 3 - (H 6 - H S ) - 211C7S1 


C7 


0 if (H 5 - H 6 )S! i 0 

1 if (H 5 - H 6 )Si < 0 
The total length of the trajectory is therefore 

TR - Ri|TRi| + R 2 |tR 2 | + «3 1 TK , 1 


(Ai2) 


(AH) 
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APPENDIX B 


FORTRAN SUBROUTINE FOR CAPTURE ALGORITHM 


SUBROUTINE NEWPST IX1,Y1*HI,R1 ,X4*Y4,R2 ,H4.X2 *Y2,H2, 

1X3* Y3*TRl*TR2*TCl»T 02*0) 

INPUTS— INITIAL POSITION* HEADING* AND TURN RADIUS Xl«Yl,n*,Kl 
FtNAL POSITION* HEADING* AND TURN RAHTUS X4*Y-*,n4»*2 
IF FINAL HEADING I SNT SPECIFIED SCT P2 TO 0. 

HI AND HA MUST BE IN RANGE C* -PI TO P? 

TWOPI«6.2B3lB53I 
l CONTINUE 

CrSH4«C0$(H4) 

SINH4*$TN(M4) 

COSHMCOSCHl) 

SINHl-SINIHl) 

C DTNIN IS MINIMUM PATH LENGTH 
0TM!N«1.EU0 
KFM*0 
KFS*1 
C 

C SELECT PATTERN TYPE 
12 CONTINUE 
C 

C RIGHT-STRAIGMT-RIGHT 
Sl«i.O 
S2'l.O 
GO T 0 20 
14 CONTINUE 
C 

C LFFT-STR AIGHT-P IGHT 

SI— 1. 

$2«l. 

K C S*2 
GO TO 20 
lfc CONTINUE 
C. 

C P IGHT-STRA IGHT-LEFT 
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d 


mmtm ittiflfa 


o o o o on o n non 


SIM. 

S2«-l. 

KFS «3 

an th 20 

i*j contthuf 

r 

C LF e T-ST°. AIOHT-l EFT 
SI — 1 . 

S?»-U 

KFS»* 

20 CrNTIMUF 


FTNP CENTERS OF TURN'S 
SI0NPl«Sl*»l 
SIGNP2«S2*R2 
XCl»Xl-Sir»NPl PSIMH1 
YCl*Yl*STGNRl*COSM 
XC2»X4-SIGNR2*SINH4 
YC2«Y4*SIGNP?*C0SH4 

0S0 IS SQUARf nF OISTANCF BETWFFN f*FMT r B s 
XN|«XC2-XCl 
VNl*YC 2-YC1 
CS0*XN 1**2*YN 1**2 

DSC IS SQUARE HP TANGENT PATH 
nSG«ST GNRI-SIGNR2 
DS0«0SC“DSG**2 

IF nsxo NO SOLUTICN 
IFfOSO.lT.O.I on 70 36 
0«SORT(PSOI 

H 2 is hfaoing of takofnt vector 
COSH 2« {n*XNi-CSG*VNl» /oso 
SINH2* (0$G*XNl*0*YNll /0S0 
H2-ATAN2! SINH2.COSH2) 

C COMPUTE LENGTHS OF TURNS 
TRI-H2-HI 

IFCPI.FO.O.I TR1*0. 

IFUTPl*Sll.GF.O. I GO TO 24 
TRI-TR l*THOPI*SI 
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24 CONTINUE 
T72=H4-H2 

I c t R? • EQ.O. I TR2*0. 

IF( CTP2*S2) .OF.O. 1 CO T n TO 
TP2=Tt?24-TWOPT*S2 
TO CONTINUE 

TD1=R l*ABS(TPl) 

T02=R2*APS(TR2I 

C OTOT TS TOTAL LENGTF CP ° ATH 
OTnT=Tni+TD2*D 
IF(KP‘«.FQ.7) GO TC 500 

FINDS IU° TO I FOUR SOLUTIONS ANO USES SH0OT=ST 


fCMOARF OTOT WITH PREVIOUS MINIMUM 
IFIOT'MN.LE.OTOT) GC TO 36 
0T M IN=OTOT 

kfm=kps 

36 CONTINUE 

CO TO (l4fl6»l8*4C)tKFS 
40 CONTI NUF 
KFS=K C M 
KPM=7 

GO TO (12»14»16»500)»KFS 
500 CONTTNUF 

COMPUTE POSITIONS OF TANGENT POINTS 
X2=XC 1+SIGNP1*SINH2 
Y?=YCl-SIGNRl*CnSH? 
X3 s XC2+SIGNP2*SINF2 
" , = YC'»-$IGNR2*C0SH2 
RFTIJPN 
END 
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Figure 2.- Calculation of capture paths 




















(a) SEQUENCE OF CAPTURES OF WAYPOINT 6 FROM DOWNWIND APPROACH. TURN 
RADil 5000 ft, INITIAL SPEED 140 knots, SPEED AT WAYPOINT 6 73 knots. 
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(b) TWO CAPTURES OF WAYPOINT 4. TURN RADII 6300 ft AND 3600 ft. INITIAL 
SPEED 100 knots, SPEED AT WAYPO»NT 4 140 knots. WIND - 20 knots. 
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SECOND TURN 


Figure 5.- Right-left-right pattern with corresponding right-straight-right 

case. 



